Properties of a diagonal 2-orbital ladder model of the Fe-pnictide superconductors 
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We study a diagonal 2-orbital ladder model ol the Fe based superconductors using the den- 
sity matrix renormalization group method. At half filling, we find a close competition between a 
"spin-striped" state and a non-collinear "spin-checkerboard" state, as well as significant nematic 
correlations. Upon finite hole or electron doping, the dominant pairing correlations are found to 
have Ai :9 (S— wave) symmetry. 



The recent discovery of iron pnictide superconductors^ 
has added to the list of materials for which the supercon- 
ducting pairing mechanism appears to be of electronic 
origin. From a theoretical point of view, these materi- 
als provide us with a rare opportunity to test our un- 
derstanding of unconventional superconducting mecha- 
nisms. Numerous models have been proposed for these 
materials. In order for these models to be more tractable, 
most authors have taken either a wealP^ or a strong^E^ 
coupling starting point (i.e. assuming that the inter- 
action strength is either much smaller or much larger 
than the bandwidth). However, there is evidence that 
the actual materials are in the intermediate coupling 
regime^K^, which is also the most difficult to treat an- 
alyticall y. In th is regime, one has to resort to numerical 
methodd 15 * 16 * 18 * ^, which are currently limited to either 
small clusters or to one dimensional systems. Despite 
these limitations, one may still hope that the important 
ordering tendencies of the system (which are presumably 
driven by short range microscopic energetics) may be ap- 
parent already for small size systems. 

In this Paper, we study a 2-orbital diagonal ladder 
model of the pnictidcs. The model is solved using the 
density matrix renormalization group (DMRG)P tech- 
nique, which enables us to study the sensitivity of the 
results to the system size in one direction. The diagonal 
geometry has the advantage that it preserves the sym- 
metry between the x and y directions, thus enabling us 
to address some of the outstanding questions of the field, 
such as the question of the gap symmetry and of nematic 
ordering. For example, the reflection symmetry of the 
model makes the distinction between Ai 9 -likc (S— wave) 
and B2g-like (D x 2_ y 2 — wave) precise. 

Model- The geometry of the model is shown in Fig. [T] 
Starting from the 2D square Fe lattice, we cut out four 
parallel chains directed along (1,-1). We then impose 
periodic boundary conditions in the transverse direction, 
such that Fe atoms separated by AR =2a (1, 1) are iden- 
tified, (a is the Fe-Fe spacing.) For each Fe site we keep 
one d xz and one d yz orbital. The Hamiltonian is written 
as H = H a + H int , where^ 
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FIG. 1: (Color online.) (a) The diagonal ladder geometry. 
Each site has two orbitals, d xz and d yz . Sites separated by 
(2a, 2a) (indicated by the arrow in the figure) are identified, 
(b) the Brillouin zone of the two-band modep'in the unfolded 
(one Fe/unit cell) scheme. The Fermi surfaces are shown, 
along with their orbital content. The diagonal lines show the 
k-states which can be accessed in the diagonal ladder. 
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where r\ = (a, a), and 
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■ ;,! — x,y creates an electron at site R 
with spin a in the d xz and d yz orbital, respectively. We 
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J2<t<t' ^L^^au' where r are Pauli matrices. In Eq. |1^ 
we have used the tight binding parameters which were 
obtained in Ref. 4 : h = -1, t 2 = 1.3, t 3 = t A = -0.85. 
(We will henceforth measure all energies in units of 
The interaction parameters in Eq. ^ were chosen to 
satisfy the constraints U — V = j J, J' = J/2 which 
follow from assuming a rotationally invariant Fe atom. 
(To check the sensitivity of the results to parameters, 
we repeated the calculations using U — V = J. The 
results did not change qualitatively.) In order to reduce 
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FIG. 2: (Color online.) (a,b) Spin patterns in 4\/2a X \/2« 
half filled (two electrons per site) systems. The arrows are 
proportional to the local (S). The diagonal dashed lines show 
the actual boundaries of the system. Zeeman fields of strength 
h — 0.5 where applied on the (0, 0) and (f , 0) sites, in the di- 
rections indicated by the bold arrows, (c) A schematic phase 
diagram at half filling. 



the number of parameters, we fixed J = \. Most of the 
calculations described below were done with [7 = 4 — 8, 
which is in the "intermediate coupling" regime (such that 
U is smaller than the overall bandwidth, but larger than 
the Fermi energy of the electron and hole pockets). 

In k-space, the diagonal ladder geometry can be 
thought of as cutting through the 2D (one Fe/unit cell) 
Brillouin zone along the lines k = k% (1, —1) + &2 (0, 1), 
where k\ € [— ^ , ^] and k% = or ^. The resulting al- 
lowed points in k-space are shown in Fig. [T|d, along with 
the Fermi surface of the 2-band model at half filling (one 
electron per orbital). The 2-orbital model has several 
well-known shortcomings: the a 2 pocket is centered at 
( - , - ) , while density functional theory calculations show 
it at (0,0), and the d xy contribution which appears on 
parts of the Fermi surface is missing. However, as we will 
discuss, the important interplay between the d xz and d yz 
is taken into account. 

The diagonal ladder has translational symmetry and 
reflection symmetry with respect to two mirror planes 
which arc formed by the z axis and a line in the (1,±1) 
directions that passes through a site. (Note that these re- 
flection operations interchange the d xz and d yz orbitals.) 
Therefore, the diagonal ladder can support a "nematic" 
phase in which the symmetry between the x and y di- 
rections is spontaneously broken. In addition, there is a 
sharp symmetry distinction between an Ai g and a Bi g 
superconducting order parameter: Ai g is even under re- 
flection through the mirror planes, and Bi g is odd. Ai g 
and B 2g (D xy ) are not distinct, since they are both even 
under reflection. However, one can still distinguish be- 
tween Ai g - and B2 g -like order parameters according to 
the relative sign of the order parameter on (1,1) and 
(1,-1) oriented bonds. 

Magnetic and nematic correlations- We begin from the 



half filled case (one electron per orbital). In order to 
study the ordering tendencies of the system, we apply 
various types of symmetry breaking perturbations at the 
edge and study how they propagate into the bulk. 

In the DMRG calculations described below, we have 
kept up to in = 6000 states in situations where both the 
number of particles and the z component of the total spin 
are conserved, and up to 3600 states in cases where one 
of these conservation laws is not present. The maximum 
truncation error was less than 4 • 10~ 4 in all cases. 

Fig. [2] shows the expectation value of the total mag- 
netization m = 5Z a (S a ) as a function of position in a 
4i/2a x \[2a system with U = 8. The total number of 
sites is 16 (32 orbitals). In these calculations, a Zeeman 
field term of the form — h R S R was applied to two sites 
near the upper left edge, R = (0, 0) and (a, 0). In Fig. [2^,, 

" (0,0) = ~h% while in Fig. 
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hz. The magnitude of 
c fields was h — 0.5. As can be seen in the figure, these 
edge fields pin very different ordering patterns. The pat- 
tern in Fig. [SJd is the "spin striped" pattern with momen- 
tum (0, f), while in 2a we find a "checkerboard" phase in 
which the magnetization on each of the two sublattices is 
orthogonal to the other, and both (— ,QJ and (0, -) mo- 
mentum components are present. This phase has been 
found in unrestricted Hartree-Fock calculations^!, and 
has been termed "orthomagnetic" (OM). The ground 
state energies in these two ordered states are equal within 
our numerical accuracy. In order to determine which 
of these states is the true groundstate, we have under- 
taken two independent approaches. We have computed 
the ground-state energy of systems with A = 8 — 16 
sites, in the presence of a bulk ordering field of strength 
h = 0.5 on each site, with a pattern of orientations which 
forces the spin-orders in Figs. [2ja,b). For U < 18, we 
have found for all system sizes studied that the energy of 
the OM pattern is lower than that of the striped pattern 
by a small amount of about 10~ 3 per site. Conversely, 
for larger U we find that the stripe pattern has a lower 
energy. We have obtained similar results from an anal- 
ysis of the "spin-nematic" order, A/r = (3(S' R ) 2 — S R ) 
measured on sites R on one sublattice with a staggered 
Zeeman field, —hS z applied to the other sublattice; pos- 
itive values of A/r arc indicative of stripe and negative 
values of OM ordering tendencies. 

The fact that the OM and the striped states are nearly 
degenerate can be understood as a consequence of the 
magnetic frustration of the ex change interactions be- 
tween the two sublattice d 9 * 10 * 13 *. In the strong coupling 
limit, U 3> 1, the 2-band model maps onto a spin- 
1 Heisenberg model with a nearest-neighbor exchange 
interaction J± and a next-nearest neighbor J 2 - In the 
regime J 2 > 0.5Ji, the classical (S — > 00) ground state 
consists of antiferromagnetically aligned A and B sublat- 
tices, while their relative orientation is completely free. 
However, 1/5* quantum fluctuations favor the striped 
state over the OM state. This accounts for the weak 
preference for stripe order for U > 18. However, for 
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FIG. 3: (Color online.) Measurements of the nematic order 
parameters [see Eq. as a function of distance from the 
origin, in &\/2a x y2a (16 site) systems. (The geometry is 
the same as in the clusters shown in Fig. [5]) (a) (7 = 4, 
(b) U = 8. In these calculations, the hopping strengths of 
the bond from (0, 0) to (a, 0) (at the upper left corner of the 
system) were 50% stronger than in the bulk. 
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FIG. 4: (Color online.) (a) Bond pairing amplitudes for a 
6y/2a x y/2a (24 site) system with U — 4, doped with 6 holes. 
Left: d xz — d xz pairing, right: d yz — d yz pairing. An external 
pairing field of strength A = 0.5 is applied to the upper left 
diagonal d xy — d xy bond [see Eq. Q], highlighted in the 
figure. The thickness of the bonds represents the amplitude 
of induced pair field on that bond. Solid (dashed) lines stand 
for positive (negative) amplitudes, respectively, (b) Detailed 
pair structure near the middle of the system. Upper left: 
d xz — d xz amplitudes, upper right: d y 
d xz — d yz amplitudes. 



somewhat smaller U , biquadratic spin interaction terms 
of the form K (R u . . . ,R 4 ) (S Rl • S R J (Sr 3 • S R J, are 
generated. Since their magnitude is of order ^ ~ J 2 /U, 
they are negligible in the large-t/ limit, but (at least for 
the parameters we have explored) they produce a weak 
preference for the OM state for U < 18. 

In the parent FeAs compounds, the ground state has a 
striped spin pattern, in contrast to our model in the re- 
alistic intermediate-coupling regime. Since the OM and 
the striped states are extremely close in energy, it is easy 
to understand how small perturbations, e.g. slightly dif- 
ferent model parameters or the coupling to the lattice, 
can stabilize the striped state relative to the OM state. 

The striped state breaks the C4 symmetry of the square 
2D lattice down to C2, and thus this state has a "ne- 
matic" component. In the diagonal ladder geometry, 



the nematic character appears as a breaking of reflection 
symmetry about the (1, ±1) directions. It is therefore in- 
teresting to calculate the nematic response of the system. 
We define the nematic order parameters: 

Afe (R) = y^([S a ,R • S ai R +tt x — S ai R ■ S a! R_ a ^]), 

a 

JV K (R) = 2X)[K<^Rda,a,B +0 *>|-Kd t a,«r > Rda,«r,R-«* 
a, a 

/Vo(R) = (Kr-t^r]). (3) 

These order parameters were measured in a calculation 
in which the hopping strength on the bond from (0, 0) 
to (a, 0) was enhanced by 50% relative to the bulk, thus 
breaking reflection symmetry about (1,-1) locally. Fig. 
[3]shows the order parameters as a function of the distance 
from the boundary at which the perturbation is applied. 

After an initial decay, Ng and Njc saturate and re- 
main nearly constant. For {7 = 4, there is some degree 
of "orbital order" No, which seems to fluctuate around 
zero. For [7 = 8, No is very small, as can be expected 
from the fact that strong repulsive interactions suppress 
both intra- and inter-orbital density fluctuations. The 
substantial nematic correlations reflect the closeness in 
energy of the spin-striped state to the ground state. 

Pairing correlations- In order to study the pairing re- 
sponse of the system, we have added a boundary pairing 
potential of the form 
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where Ri = (0,0) and R2 = (a, a). This term can be 
thought of as a proximity coupling to a bulk supercon- 
ductor. Note that since the pairing term is applied only 
to the d xz orbital, it couples to any (singlet) supercon- 
ducting order parameter. 

Fig. [4^i shows the induced bond pair amplitudes 

0R,R' = (<4,t,R d Ia,R' - 'C. li'C-.li/ + H - C -> for a 

6\/2a x y/2a (24 site) system with A = 0.5 and U — 4. 
With the pair field term Q, the number of electrons in 
the system is not conserved. Here, we choose the chemi- 
cal potential \x such that the average number of electrons 
in the system is close to 42, which corresponds to a hole 
doping of n = 0.25 per Fe site. The pair amplitudes decay 
slowly with the distance from the edge. Near the mid- 
dle of the system, the pair structure is shown in Fig. [4Jd. 
One can see that under reflection about the (1, ±1) direc- 
tions, the pair wave function is even, indicating A lg -like 
pairing. The Ai s structure appears already in smaller 
systems (down to 2^/2a x ^/2a). The same pairing oc- 
curs in electron-doped systems and even in the undoped 
system, although here the pairing response is consider- 
ably weaker. For U = 8, the pairing symmetry is still 
Ai g , although the pairing amplitude decays considerably 
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FIG. 5: (Color online.) Induced pair amplitudes on diagonal 
bonds in three 4y / 2a x \/2a calculations, in which Ai l9 (□), 
Bi i9 (o) and triplet (❖) pair fields are applied to the edge (see 
text). Solid (dashed) lines correspond to d xz — d X z (d yz —d yz ) 
pairing, respectively. 

with various pairing potentials applied at the edge, such 
that they excite different pairing symmetries selectively. 
For Ai g , Eq. Q was used. For Bi g , we have added an 
additional term for the d yz orbital, with an opposite sign. 
This term does not couple to the A l£( pairing symmetry 
at all, because it is odd under reflection about (1,-1). 
In addition, we have applied a triplet pairing term [in 
which the sign in the square brackets of Eq. Q is re- 
versed]. Strictly speaking, the B2 9 (D xy — like) symmetry 
is not distinct from A± g . However, one can still think 
about a "E^g— like" pair structure in which the sign of the 
pair amplitude on (1,1)— and (1, —1) —oriented bonds 
is opposite. To couple mostly to B2 ff , we have added 
to Eq. Q a bond pairing term on the (0,0) — (a, —a) 
bond with an equal magnitude and opposite sign. Alter- 
natively, we have applied a singlet d xz — d yz pair field on 
the (0, 0) — (a, 0) bond, which also favors B2 S -like pairing. 
Comparing the response of the different pairing symme- 
tries, the Ai g pairing clearly decays much more slowly 
than the triplet and the B ls pairings. In the B 2g cal- 
culations (not shown), the induced pairing structure is 



B2 9 — like very close to the edge, but changes its charac- 
ter to Ai g — like further away into the bulk. 

Discussion- In conclusion, the two-orbital diagonal 
ladder model shows several robust features, which ap- 
pear already at very small system sizes. In electron and 
hole doped systems, there is a clear tendency to form 
Ai g (S— wave) superconducting order. All other forms 
of order are very weak. For the undoped system, a 
strong tendency toward antiferromagnetic ordering is ob- 
served. However, two forms of magnetic order are in close 
competition with each other: unidirectional "stripe" or- 
der, of the sort found in experiment, and non-collinear 
"spin-checkerboard" (OM) order, which was found also in 
mean-field calculations^. In our model, the spin-striped 
state is stabilized for U > 18, and the OM is stabilized 
for smaller values of U. The near degeneracy of these 
two states implies that small terms can tilt the balance 
one way or the other, which raises the possibility that the 
OM state may be stabilized in some member of the Fe- 
pnictide or chalcogen families. Finally, in the undoped 
system there is a strong tendency towards nematic or- 
der. This order is associated with the expectation values 
of bond operators (e.g. the local kinetic energy and spin- 
spin correlations), rather than the difference of onsite 
orbital occupations which were found to be small. 
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